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We discuss and compare different approaches to include gluon transverse momenta for 
heavy quark-antiquark pair and meson production. The results are illustrated with the 
help of different unintegrated gluon distributions (UGDF) from the literature. We com- 
pare results obtained with on-shell and off-shell matrix elements and kinematics. The 
results are compared with recent experimental results of the CDF collaboration. 



1. Introduction 

The heavy quark-antiquark production in hadroproduction is known as one of the 
crucial tests of conventional gluon distributions within a standard factorization 
approach. At high energies one tests gluon distributions at low values of longitu- 
dinal momentum fraction. Standard collinear approach does not include transverse 
momenta of initial gluons, the method to include transverse momenta is k t - factor- 
ization approacrPEEl In the first step of our approach the single particle spectra of 
charmed quarks and antiquarks are obtained assuming gluon-gluon fusion . Differ- 
ent unintegrated gluon distributions from the literature are used^SlZElSI To obtain 
the single particle spectra of meson from those of quarks/antiquarks a standard 
hadronization procedure with Peterson fragmentation function is applied. Conclu- 
sions about unintegrated gluon distributions are drawn. 
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Fig. 1. Leading-order diagrams for heavy quark - antiquark production. 



2. Heavy quark production 

Let us consider the reaction h\ + h% — * Q + Q + X , where Q and Q are heavy quark 
and heavy antiquark, respectively. 

In the leading-order approximation within collincar approach the triple-differential 
cross section in rapidity of Q (yi), in rapidity of Q (1/2) and transverse momentum 
of one of them (p t ) can be written as 

da 1 



\ " 



xipi(xi,n ) X2Pj{x 2 ,n ) \Mi 3 \ 2 . (1) 



dyidy2d 2 pt 16n 2 s 2 ^r! 

Above pi(xi,pr) and Pj(x2,p?) are familiar (integrated) parton distributions in 
hadron hi and I12, respectively 

The parton distributions are evaluated at: x\ — ^= (exp(j/i) + exp(y2)), £2 = 
^= (exp(— j/i) + exp(— 2/2))- The formulae for matrix element squared averaged over 

initial and summed over final spin polarizations can be found e.g. in 

If one allows for transverse momenta of the initial partons, the transverse mo- 
menta of the final Q and Q no longer cancel. Formula JIJ can be easily generalized 
if one allows for the initial parton transverse momenta. Then 

da ■sr-y f d 2 Ki t t d 2 K2,t 



dy\dy2d 2 p\jd 2 p2,t J tt tt 1Qtt 2 {x\X2s) 2 



\Mi 



2 



S 2 (k m + K 2 ,t -Pi,t ~P2,t) /i(aii, «?,t) fj( x 2,i4,t) ) ( 2 ) 

where now fi(xi,K 2 t ) and /^(a^, ftf t) are so-called unintegrated parton distribu- 
tions. The extra integration is over transverse momenta of the initial partons. The 
two extra factors 1/7T attached to the integration over d 2 K\^ and d K2,t instead 
over <ifq t and c?K2 jt as in the conventional relation between unintegrated and in- 
tegrated parton distributions. The two-dimensional delta function assures momen- 
tum conservation. Now the unintegrated parton distributions must be evaluated at: 
xi = ^-exp(yi) + ^exp(y 2 ), x 2 = ^exp(-yi) + ^exp(-y 2 )- In general, 



*Based on a talk given by M. Luszczak at MESON2006 (Krakow, June 9-13, 2006) 
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the matrix element must be calculated for initial off-shell partons. The correspond- 
ing formulae for initial gluons were calculated in 1^21 ( see a i so GO) j n the present 
paper for illustration we shall compare results obtained for both on-shell and off- 
shell matrix elements. 




Fig. 2. Inclusive dcr/dpt for charm/anticharm production at W = 1.96 TeV for different UGDF. 
The results with on-shell kinematics are shown in panel (a) and results with off-shell kinematics 
in panel (b). In this calculation both factorization and renormalization scales were fixed for 4 m\. 



In Fig. 2 we collected results for da/dp t obtained with different unintegrated 
gluon distributions from the literature. In this case consequently the off-shell ma- 
trix element and off-shell kinematics were used. The GBW gluon distribution leads 
to a much smaller cross section. The KL gluon distribution produces the hardest 
pt spectrum. Rather different slopes in transverse momentum of c (or c) are ob- 
tained for different UGDFs. This differences survive after convoluting the inclusive 
quark/antiquark spectra with fragmentation functions. Thus, in principle, precise 
distribution in transverse momentum of charmed mesons should be useful to select 
a "correct" model of UGDF. A detailed comparision with the experimental data 
requires, however, a detailed knowledge of fragmentation functions. 

The inclusive spectra are not the best observables to test UGDF H2. Let us 
come now to correlations between charm quark and charm antiquark. 

In Fig. 3 we compare results for different unintegrated gluon distribution from 
the literature. Quite different results are obtained for different UGDFs. The non- 
perturbative GBW glue leads to strong azimuthal correlations between c and c. 
In contrast, BFKL dynamics leads to strong decorrelations of azimuthal angles of 
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Fig. 3. Azimuthal angle correlations for different UGDFs in the literature: Kwiecifiski (solid), 
BFKL (dash-dotted), GBW (dashed) and KL (dotted). 



charm and anticharm quarks. The saturation idea inspired KL distribution leads 
to an local enhancement for (f> c5 which is probably due to simplifications made in 
parametrizing the UGDF. In the last case there is a sizeable difference between the 
result obtained with on-shell (left panel) and off-shell (right panel) matrix elements. 
All this is due to an interplay of the matrix element and the unintegrated gluon 
distributions. 



3. From unintegrated parton distributions to meson production 

The inclusive distributions of hadrons are obtained through a convolution of inclu- 
sive distributions of heavy quarks/antiquarks and Q — > h fragmentation functions 

dcr(y h ,pt, h ) _ f 1 dz D da£ g ^ Q (y Q ,pt,Q) 
dy h d 2 p t , h J z 2 Q ^ h dy Q d 2 p tl Q VQ=Vh 

Pt,Q=1>t,h / * 

In the present paper we have used so-called Peterson fragmentation functions with 
parameters from last issue of PDG (Particle Data Group). We have neglected a 
possibility that charmed mesons are produced from light quarks and/or gluons. 
This approximation should be better for heavier quarks/mesons. 

In the present analysis we show dependence of the total cross section on trans- 
verse momenta for D* + production 

dcrjpt) = I" 1 dy da(y,p t ) _ ^ dajy = 0,p t ) 
dp t J -i dydpt dydpt 
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Fig. 4. Inclusive da/dpt for D*+ production at W = 1.96 TeV for different UGDF. The results 
with ce s = a s (4m^) is in the left panel and the results with a s = a B (K^ t ) or a s (K^ t ) in the right 
panel. 



In Fig. 4 the theoretical results are compared with experimental data of the CDF 
collaboration El. "We collected results obtained with different unintegrated gluon 
distributions from the literature. 

We have made calculations for two different choices of renormalization scale. 
The results obtained with a s {Am 2 c ) (left panel) are below the experimental data, 
while results with a s = a s (n\ t ) or a s («2t) (see^L*") better describe the data. There 
is also a dependence on UGDF used in the calculation. Even if the second choice is 
made, there seem to be a deficit of the cross section. 

It is not clear to us if the deficit is due to the omission of q — > D* or g — > D* 
fragmentation, or due to higer order terms in heavy quark production. 

Similar situation is for B + production, as shown in Fig. 5. 

4. Summary 

Inclusive cross section for heavy quark - antiquark and heavy mesons in proton 
-(anti)proton collisions have been calculated in the /c t -factorization approach. We 
have compared quantitatively different methods to include gluon transverse mo- 
menta and their effect on inclusive spectra as well as on cc correlations. Different 
UGDF from the literature were used. 

The inclusive spectra of heavy mesons were obtained via convolution of heavy 
quark spectra with so-called Peterson fragmentation functions. The results depend 
on the choice of the factorization scale. The "best" results are obtained with the 



February 2, 2008 6:56 WSPC/Guidelines-IJMPA 



luszczak 



6 M. Luszczak and A. Szczurek 




Fig. 5. Inclusive dcr/dvt for _B+ production at W = 1800 GeV for different UGDF. The experi- 
mental data are from RefEU The results with a s = a s (4m^) is in the left panel and the results 
with a a = a a {K^ t ) or a a {K^ t ) in the right panel. 



BFKL UGDF. We observe a deficit of the cross section for almost all UGDFs used. 
It is not clear to us if the deficit is due to neglecting some hadronization components 
or due to next-to-leading order effects. 
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